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Influence of Wheel Load Shape on Vertical Stress Reaching 
Subgrade through an Aggregate Layer 


KAREN S. HENRY 


INTRODUCTION 

The U.S. Army design procedure to stabilize low- 
bearing-capacity soil with geotextiles for traffic in¬ 
volves placing aggregate on geotextiles (U.S. Army 
and U.S. Air Force 1995). It is based on the assump¬ 
tion that the applied surface load (the wheel load) is in 
the shape of a circle (e.g., Henry 1999). The vertical 
stress that reaches the subgrade from the applied wheel 
load is then estimated on the basis of the assumption 
that the aggregate is an elastic, homogeneous, isotro¬ 
pic half-space and, therefore, that the stress distribu¬ 
tion can be estimated by using the Boussinesq method 
(e.g., Newmark 1942). The stress distribution deter¬ 
mined by the Boussinesq approach is assumed to be 
accurate for the prediction of stresses distributed 
through compacted crushed rock and through asphalt 
(e.g., Barenberg et al. 1975, Yoder and Witzcak 1975). 
However, it is not clear that a circular area accurately 
predicts stress for non-circular wheel load shapes. For 
example, a common configuration of dual wheels on a 
single axle has been modeled as a rectangular area (e.g., 
Giroud and Noiray 1981). To quantify the extent to 
which the shape of the applied load influences the stress 
that reaches the subgrade through the aggregate layer, 
I analyzed the influence of the shape of the wheel load 
on the maximum vertical stress at depth predicted by 
the Boussinesq method. 

Since it is the maximum vertical stress (i.e., the 
stress beneath the center of the loaded area) that is used 
in the design procedure, I examined the influence of 
the shape of the applied load on the maximum vertical 
stress predicted by the Boussinesq approach. If sig¬ 
nificant differences in the maximum vertical stresses 
below uniformly loaded rectangles and circles were 
found, then the wheel load shape (i.e., rectangle or 
circle) should be accounted for in the design proce¬ 


dure. A comparison of the maximum vertical stress 
beneath circular and rectangular loads as a function of 
depth indicated that the differences were insufficient 
to warrant a change in the design procedure for rect¬ 
angular load geometries that have length/width ratios 
of 3 or less, and, furthermore, that assuming that the 
load is circular is conservative. 

The details of the analysis are presented fully in 
this report. The techniques may be used to examine 
other load geometries, including those of tracked ve¬ 
hicles. 


STUDY OBJECTIVES 

The goal of this study was to determine whether 
the shape of the wheel load at the ground surface sig¬ 
nificantly influences the maximum vertical stress at 
depth (i.e., that reaches the subgrade) given that the 
aggregate layer behaves as a linearly elastic material. 
If a significant difference were to be found between 
the vertical stresses at depth applied by uniformly 
loaded circles versus rectangles, then appropriate 
changes would be made in design guidance developed 
for vehicles that apply loads that are more accurately 
modeled as rectangles than circles. 

I used the Boussinesq equations to estimate maxi¬ 
mum vertical stresses for uniformly loaded circles 
(Newmark 1942) and rectangles (Newmark 1935). 1 
assumed unit loads and areas, and considered a large 
range in length-to-width ratios for rectangular wheel 
loads so that the range of wheel load shapes applied 
by military vehicles was well represented. The longest 
rectangle that I studied, at L — 6B, where L is the con¬ 
tact length and B is the contact width, is too extreme to 
model wheeled vehicles. The range of normalized 



depths studied represents a range of typical aggregate 
thicknesses above geotextiles. The technique used can 
easily be applied to other length/width ratios and other 
normalized depths. 


VERTICAL BOUSSINESQ STRESS 
BELOW A RECTANGLE 


The vertical stress at depth below the comer of a 
uniformly loaded rectangle of dimensions L x B is es¬ 
timated as (Fig. 1; Newmark 1935) 

0 = jny 2 A o V o V o +1 (la) 

z ~ 47r f 0 2 + a 2 v 2 


tan -1 - 0 0 for ^ > Z 

k 2 -a 2 

o o y 


„ _ _w_ 2AK KT ± 1 I 
r/H 1 ' C 


tan 1 — + 71 for A a < Z 

Vo 2 -A 2 ) 


where 

w = applied uniform pressure 
Z = depth below the corner of the rectangle 
for which the vertical stress is desired 



Figure 1. Notation used in the estima¬ 
tion of vertical stress at depth below 
the corner and center of a uniformly 
loaded rectangle. 

and 

V 2 _ X 2 + Y 2 + Z 2 
Z 2 

Thus 1 , or the ratio of the vertical stress below the 

z 5 

center of a rectangle at depth Z to the stress applied on 

the rectangle (L x B) at the surface, is equal to 


2 _ I 2 +B 2 +Z 2 
° ~ Z 2 


Note that V Q is the ratio of the radius vector from point 
Z to the corner of the rectangle to the depth of point Z 
To obtain the stress below the center of the rectangle 
at depth Z, I used superposition and divided the rect¬ 
angle into four smaller rectangles, such that 
„ L 


andZ and V are defined as follows: 
a-XY 


*This form of the equation is not published by Newmark 
(1935), but it is required to avoid predicting tension (e.g., 
Bowles 1988). 
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VERTICAL BOUSSINESQ STRESS 
BELOW A CIRCLE 

The vertical stress at depth Z below the center of a 
uniformly loaded circle of radius r, divided by the ap¬ 
plied stress, is estimated as (Fig. 2; Newmark 1942) 


2 




y 



Figure 2. Uniformly loaded circle, 
showing notation used in the esti¬ 
mation of vertical stress at depth 
below the center of the circle. 
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I used eq 2 and 3 for a unit uniform load and a unit 
area to compare the stress reduction at various nor¬ 
malized depths estimated by each equation. 1 varied 
the length/width ratio of the rectangles from L = B to L 
= 6 B and the normalized depth Z from 0.05 to 2.0. The 
range in these parameters reasonably represents the 
geometry of loads applied by heavy military vehicles. 
An example of estimating the length-to-width ratio of 
a wheel load and the values of V and r for a military 
vehicle follows. 


WHEEL LOAD GEOMETRY 
OF THE M1000 TRAILER 

On-highway trucks typically have L.B ratios of 1.4 
and off-highway trucks have L:B ratios of 2.0 (Giroud 
and Noiray 1981). The Ml000 trailer of the Heavy 
Equipment Transporter System (HETS) (Fig. 3) likely 
has the highest length/width ratio of a wheel load ap¬ 
plied by a U.S. military vehicle. Each axle line of this 
trailer has eight wheels (two sets of dual wheels on 
each end) (Fig. 4). If we assume that the soil between 
the tires aligned on one side of the axle will not fail 
(e.g., Giroud and Noiray 1981), then the load carried 
on that side of the axle can be assumed to be distrib¬ 


uted over a rectangle (Fig. 4). Based on the informa¬ 
tion given in Figure 4, the L:B ratio of the uniformly 
loaded rectangle is approximately 3.9. Note that for a 
1-m depth of aggregate, V is equal to 1.16 for this ve¬ 
hicle. 


ANALYTICAL STUDY AND RESULTS 

I programmed eq 2 and 3 into a Mathcad 8 pro¬ 
grammable worksheet (Mathsoft 1998) and set the area 
of pressure application equal to one. The results gen¬ 
erated include the fraction of the applied vertical stress 
below the rectangle and the circle, 7 ect and / irc , respec¬ 
tively, and the ratio I 11. . A factor F represents the 
L.B ratio that was easily changed to examine its influ¬ 
ence on results. Appendix A contains a hard copy of 
the Mathcad worksheet used, tabulated results, and a 
copy of the Mathcad worksheet for the reader’s use. 

The ratios of the stress reduction beneath the center 
of the rectangles to the stress reduction beneath the 
center of a circle are shown in Figure 5. The maxi¬ 
mum vertical stress beneath the circle is always greater 
than that beneath the rectangles, including the special 
case of a square. A circle minimizes the ratio of perim¬ 
eter to area; thus, the circle maximizes the influence 
of surface pressure at depths below the center. There¬ 
fore, applying the assumption that the surface load is 
circular is conservative compared to assuming other 
shapes. Furthermore, for length/width ratios of 3 or 
less, the maximum vertical stress applied by a rect¬ 
angle is at least 85% of that applied by a circle. The 
minimum ratio of the stress beneath a rectangle to a 
circle for the conditions analyzed occurs at a normal¬ 
ized depth of about 0.6, and is about 0.65 for a length/ 
width ratio of 6. As mentioned earlier, the length/width 
ratio of 6 is outside the range of load shapes for wheeled 
U.S. military vehicles. 



Figure 3. M1000 semi-trailer, downloaded from 
USATACOM (1999). 
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A = 178 mm (7 in.) 

B = 305 mm (12 in.) 

C = 305 mm (12 in.) 

D = 85.7 mm (3.375 in.) 


Figure 4. Dimensions of wheels on one end of an 
axle line of the Ml 000 trailer.* 

DISCUSSION 

Based on the results of the analytical study, there is 
no compelling reason to assume any wheel load shape 
other than circular. However, if a vehicle does produce 
wheel loads than can be approximated as uniformly 
loaded rectangles with length/width ratios of about 4 
or greater, this question should be reconsidered for pur¬ 
poses of modifying design guidance that may be too 
conservative. In this case, the spreadsheet shown in 
Appendix A can be used to calculate the stress reduc¬ 
tion at depth below the rectangle simply by changing 
the F value that is input at the top of the spreadsheet. 
Keep in mind that this approach assumes that the 
Boussinesq stress distribution provides a reasonable 
estimate of vertical stress at depth (i.e., the aggregate 
behaves in a linearly elastic manner and the stress is 
uniformly applied). 

Although not the topic of this study, whether the 
aggregate that is placed on top of the geotextile behaves 
as a linearly elastic material is important; related to this 
is whether the vertical stresses are accurately predicted 
by the Boussinesq method. The aggregate type and qual¬ 
ity significantly influence the stress distribution at depth. 
McMahon and Yoder (1960) found that pressures in 
homogeneous, compacted clay soil subgrades followed 
the Boussinesq pattern, but were up to 18% higher than 
those predicted by Boussinesq at shallow depths (e.g., 
less than the diameter of the loaded area). Further, the 
presence of compacted crushed limestone base above 
the clay subgrade reduced the vertical stresses in the 
subgrade immediately below the interface to values less 
than those predicted by Boussinesq. Other research in¬ 
dicated that there is a stress concentration in soil layers 


*Personal communication with R. Dunbar, Test Director, 
Yuma Proving Grounds, Yuma, Arizona, 1999. 


containing significant quantities of sand compared to 
other aggregates. For example, a 45-kN (10-kip) load 
was applied at 690 kPa (100 psi) to a subgrade through 
various aggregates, with the result that the stress reach¬ 
ing the subgrade through a 0.6-m (24-in.) layer of sand 
was about twice that of the stress through a layer of 
crushed limestone (Hemer 1955). Sowers and Vesic 
(1962) found that vertical stresses in a 203-mm- (8- 
in.-) thick base consisting of 95% subangular quartz 
sand and 5% asphalt, overlain by a 152-mm- (6-in.-) 
thick asphalt pavement, were up to 70% higher than 
the vertical stress predicted by Boussinesq and higher 
than that in other types of base—many of which had 
lower vertical stresses than those predicted by 
Boussinesq. Thus, more effort should focus on accu¬ 
rately modeling stress distribution through the vari¬ 
ous aggregates that the military may use, indeed may 
be forced to use, in constructing geotextile-reinforced 
unsurfaced roads. 


CONCLUSION 

For all wheeled U.S. military vehicles, the wheel 
load may safely and accurately be assumed to be cir¬ 
cular if the aggregate that is loaded can be modeled as 
a linearly elastic material (e.g., compacted, crushed 
rock). Because the U.S. Army design procedure to use 
geotextiles overlain with high-quality aggregate to sta¬ 
bilize low-bearing-capacity soil for trafficking with 
wheeled vehicles (U.S. Army and U.S. Air Force 1995) 
assumes that the wheel load applied at the surface is 
in the shape of a circle, no changes are required in this 
approach to modeling wheel loads. However, if a ve¬ 
hicle does produce loads than can be approximated as 
uniformly loaded rectangles with length/width ratios 
of about 4 or greater, (e.g., tracked vehicles) this ques- 



Figure 5. Ratio of the reduction in applied stress 
with depth for a unit rectangle to a unit circle for 
various length/width ( L\B ) ratios of the rectangle. 
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tion should be reconsidered in order to save aggregate 
material. Furthermore, the Mathcad spreadsheet as 
shown in Appendix A can be used to calculate the stress 
reduction at depth below the rectangular loads. 

More important than wheel load shape is whether 
the aggregate that is placed on top of the geotextile 
can be accurately modeled as a linearly elastic mate¬ 
rial, and, therefore, whether the vertical stresses are 
reasonably approximated by the Boussinesq method 
(as the current design method assumes). Indeed, stress 
concentrations in sand and in clay soil compared to 
Boussinesq predictions have been noted. Thus, if low 
quality aggregate composed of rounded gravel, con¬ 
taining sand or fines, or both, must be used in geotextile 
reinforcement of low-bearing-capacity soil, this ques¬ 
tion should be addressed. 
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APPENDIX A: ANALYTICAL STUDY OF MAXIMUM VERTICAL STRESSES 
UNDER UNIFORMLY LOADED RECTANGLES AND CIRCLES 

This worksheet,nrmlshapestudymatrix.mcd, prepared by Karen S. Henry of The U S. 
Army Cold Regions Research and Engineering Laboratory, Hanover, NH, compares the 
distributions of vertical stress at depth of uniformly loaded circles and rectangles by the 
respective Boussinesq equations. F can be varied to study the influence of the 
length/width ratio, z may also be varied to determine the stress at various depths. 


Assume unit Areas: I Area 

Area := 1 r:= I- r = 0.564 F := 6 

yj % 

F is the ratio of the length, L to the width, B of the rectangle 
unit area. Note that z is a vector. 

i:= 0.. 19 

B(F) := L(F) := F-B(F) z. := + 0.1 L(F)-B(F) = 1 


Divide the rectangle into four sections: 


X(F) := 


UF) 


Y(F) := 


B(F) 


X(F) = 1.225 


2 2 
Define V and A for use in Newmark's rectangle equation. 


Y(F) = 0.204 


V(z,F) 


2 2 2 
X(F) + Y(F) + z 


\1 z 2 


A(z,F) := 


X(F)-Y(F) 

2 

z 


This expression, of the fraction of the applied load at z, is for the 
sum of four rectangles of dimension X by Y. 


Irectl(z,F) 


Irect2(z, F) 


f 

2-A(z,F)-V(z,F) 

2 ^ 
V(z,F) + 1 

+ atari | 

f 2-A(z,F)-V(z,F) > 

, V(z,F) 2 + A(z,F) 2 

V(z,F) 2 j 

v V(z,F) 2 - A(z,F) 2 J_ 


1 

Y 2 ) 

2-A(z, F)-V(z,F) V(z,F) + 1 

+ atanj 

f 2-A(z,F)-V(z,F) > 

+ 7C 




n 

_,V(z,F) 2 +A(z,F) 2 V(z,F) 2 j 


v V(z,F) 2 -A(z,F) 2 v 



Irect(z,F) := if(A(z,F) > V(z,F),Irect2(z,F),Irectl(z,F)) 


Icircfz) := 
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Results: 



12.457 

6.288 

4.258 

3.261 

2.677 

2.298 

2.036 

1.846 

1.704 

1.594 

1.508 

1.439 

1.383 

1.337 

1.298 

1.266 


5.642 

2.821 

1.881 

1.41 

1.128 

0.94 

0.806 

0.705 

0.627 

0.564 

0.513 

0.47 

0.434 

0.403 

0.376 

0.353 


0.1 

0.2 

0.3 

0.4 

0.5 

0.6 

0.7 

0.8 

0.9 

1 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 


I recti z., F 
i , 


0.961 

0.824 

0.675 

0.556 

0.465 

0.396 

0.343 

0.299 

0.264 

0.235 

0.21 

0.189 

0.17 

0.154 

0.141 

0.128 


Icirc(z) = 


0.995 

0.963 

0.897 

0.807 

0.708 

0.613 

0.528 

0.454 

0.392 

0.339 

0.296 

0.259 

0.228 

0.202 

0.18 

0.161 


Irect|z.,Fj 


Icirc^z.j 


0.967 

0.856 

0.753 

0.689 

0.657 

0.646 

0.649 

0.659 

0.674 

0.692 

0.71 

0.729 

0.747 

0.764 

0.781 

0.796 


Note that the augment function pastes matrices next to each other. The arrays 
generated by the computations that were not already vectors have to first be "vectorized 
with the vectorize function (z and functions of z are vectors, but not functions 
containing both z and F). This is done so that all elements in the matrix generated 
by the results can be exported to a spreadsheet. 



f 

f 

f 

f 

( --^ 

v 

\\\ 

M := augment 

z, augment 

- — ^ 

V(z,F), augment 

r 

-, augment! 

- > 

Irect(z, F), augment 

T • f \ Irect(z,F) 
Icirc(z), 







V 

^ Icirqz) ) 

). 

))) 



Icirc^Zj) 
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